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Problems with the Quenched Approximation in the Chiral Limit 

Stephen R. Sharpe'" * 

^Physics Department, University of Washington, Seattle, WA 98195, USA 

In the quenched approximation, loops of the light singlet meson (the n') give rise to a type of chiral logarithm 
absent in full QCD. These logarithms are singular in the chiral limit [Q-h], throwing doubt upon the utility of 
the quenched approximation. In previous work, I summed a class of diagrams, leading to non-analytic power 
dependencies such as {^tp) oc rUq^^'^^^^'' I suggested, however, that these peculiar results could be redefined 
away. Here I give an alternative derivation of the results, based on the renormalization group, and argue that 
they cannot be redefined away. I discuss the evidence (or lack thereof) for such effects in numerical data. 



1. INTRODUCTION 

Many simulations of lattice QCD use the so- 
called "quenched" approximation, in which the 
fermion determinant is left out of the measure. 
Since this approximation is likely to be with us 
for a few years yet, it behooves us to attempt to 
understand all we can about the peculiarities of 
the quenched theory. 

I focus here on the effects of the extra light 
pseudo-Goldstone boson (PGB) present in the 
spectrum of quenched theory. I assume degen- 
erate light quarks {Nf each of mass m^) in which 
case the extra PGB is the flavor singlet state. By 
analogy with QCD, I call this the 77'. 

The 77' is a PGB in the quenched theory because 
the series of diagrams shown below are absent. 
The first of these diagrams is, however, present 

+ •) ( -) ( + . . • 

in loop graphs, and introduces a new dimension- 
ful parameter. This changes the power counting 
rules of chiral perturbation theory, and leads to 
enhanced chiral logarithms in quenched QCD. 

This effect, first pointed out in Ref. |^], 
was first calculated by Bernard and Golterman 
(BG)[|. They developed a chiral La grangian for 
quenched QCD, containing ghost states to cancel 
the effects of the quark loops. They found, e.g. 

ml/2fnng ^ [1 - 2,51n(m^VA) + 0(mj] , (1) 
S = ml/iNfWn^f^) , (2) 

where A is the ultraviolet cutoff, = 93MeV 
and /i are parameters in the chiral Lagrangian, 
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and mo is the scale characterizing the rj' two- 
point function. These parameters are discussed 
further below. In full QCD, the term propor- 
tional to S is absent, and so = 2^mq in the 
chiral limit. Chiral logarithms, proportional to 
r7iqln(mg), vanish in this limit. (These are rep- 
resented by the shorthand 0{mq) in Eq. |l|.) In 
quenched QCD, on the other hand, the logarithm 
diverges in the chiral limit. 

Using the quark-line method I obtained, inde- 
pendently, the same results for the enhanced chi- 
ral logs from rj' loops js],^. (I also noted that 
there are a number of quantities (e.g. f-„, Bk) 
which, at one-loop, and for degenerate quarks, 
are not affected by -q' loops.) In order to address 
the question of the approach to the chiral limit, 
when the term (51n(TO^') gets large, I summed the 
leading logarithms proportional to powers of this 
quantity, with the result 

ml/2txmq = (2//m,/A2)-V(i+^) + 0(m,) . (3) 

The divergence has been enhanced from a loga- 
rithm to a power law. 

This result appears to violate PCAC relations 
such as f^rn^ — —2mq{-ipip). This is not so, how- 
ever, since there is a whole tower of equations 
similar to Eq. |^, e.g. 

~m/{fif^) = (2Mm,/A2)-^/(i+^)+0(m,) .(4) 

These results are of the same form as those de- 
scribing critical behavior close to a second order 
fixed point. At first sight this is puzzling, because 
no such fixed points (with non-trivial exponents) 
are known in four dimensions. The puzzle is re- 
solved by recalling that this lore applies to critical 
theories which are unitary, and thus presumably 
not to the massless limit of quenched QCD. The 
similarity to critical behavior does, however, sug- 
gest that it should be possible to recast the results 
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Eqs. I and | into a renormalization group (RG) 
framework. 



2. RG DERIVATION 

It is simplest to use the quenched chiral La- 
grangian of Bernard and Goherman, £bg- The 
symmetry of £bg is the graded group [/(3|3)l x 
f/(3|3)i^, and there are extra pseudoscalar and 
fermionic PGBs. It turns out, however, that 
for ampUtudes involving only the physical PGBs 
as external particles, and for external momenta 
p <C 47r/^ , one need only consider part of Cbg ■ 

Cbg = Tl■[^^7^^^7^]-l^imgf^^T[U + U^ 

+ \mlTl^ + C!^^ . (5) 

The pion field here includes the ry': 

8 

TT = ^ TToTa , "Yx\TaTh\ = \5ab , U = exp[2j7r/ /^] 
0=0 

£bq includes various types of terms which can 
be neglected since they do not give rise to en- 
hanced chiral logarithms. First, there are other 
leading order terms (i.e. of 0{p^ , mq)) containing 
the additional PGB fields, which have the effect 
of canceling various diagrams involving the phys- 
ical PGB's. These cancellations are discussed in 
detail in Refs. ||,^. Most important here are 
the restrictions on the insertions of the ij'^ term. 
While this looks like a mass term, it can in fact 
only be inserted once on each internal ij' line, be- 
cause of cancellations with terms from Cqg- By 
contrast, in full QCD does contribute to the 
r]' mass, m'^y — 2^mq + TOq. A standard phe- 
nomenological analysis using the physical meson 
masses, and allowing for non-degenerate quarks, 
gives mo « 0.9GeV in full QCD. 

Using this value as an estimate in the quenched 
approximation, one finds S ~ 0.2. This esti- 
mate is unreliable for two reasons. First, pa- 
rameters in the quenched and full chiral La- 
grangians need not be the same. Second, I have 
excluded a possible wavefunction renormalization 
term ^{A — l)d^ri' d^rj' , which reduces the esti- 
mate of (5 by a factor of A. This term does not, 
however, contribute enhanced chiral logarithms, 
because it does not introduce a new scale. It sim- 
ply makes the estimate of 5 less certain. 



also includes the remainder of the full ki- 
netic term 

-^SlTr[d^Ud^U^] = Tv[d^^d^Ti] + 0(7r4) . 

Only the 0(7r^) term need be kept, because the 
higher order vertices do not contain the 77' field, 
and do not contribute enhanced chiral logarithms. 
On the other hand, the full mass term must be 
kept, since the 77' is contained in all terms. 

Finally, C'^q contains 0(p*) terms which are 
suppressed in the chiral limit. 

I now turn to the calculation. In full QCD, 
power counting shows that logarithms always 
come with a factor of the quark mass or the ex- 
ternal momenta-squared, e.g. 

ml|{2^lmq) = 1 + 2to2 in(rnJA)/(A/16^V') ■ 

The logarithmic cut-off dependence can thus be 
absorbed by the coefficients of the ©(p"*) part of 
the Lagrangian, while the leading order coeffi- 
cients are unaffected. 

If one uses a simple momentum cut-off, rather 
than, say, dimensional regularization, then there 
are also quadratic divergences. In full QCD, these 
do affect the leading order coefficients. One must 
introduce a scale dependence in the bare decay 
constant, /(A), so that physical quantities are in- 
dependent of the cut-off. These divergences are, 
however, absent in the quenched theory, so /tt re- 
mains independent of scale. 

In quenched chiral perturbation theory power 
counting is changed because of the additional 
dimensionful parameter mg. Insertions of the 
77'^ vertex into diagrams that were previously 
quadratically divergent produces logarithmically 
divergent integrals, giving factors of ^ ln(TO^/ /A). 
For example. Fig. lb leads to Eq. |l|. Since 5 does 
not vanish in the chiral limit, the scale depen- 
dence must be absorbed into the leading order 
coupling: ^ /^(A). There are also divergent 
corrections to mg, but these can be shown to be 
suppressed by powers of TO^/(47r/7r)^. 

In Ref. 1^ I summed all diagrams contain- 
ing leading logarithms, i.e. those proportional to 
[(51n(m^/)]" with no additional factors of or 
5. Two crucial points allowed this. First, in order 
to get a leading logarithm one must insert each 
ry'^ vertex in an initially quadratically divergent 
loop. Furthermore, every loop in the initial di- 
agram must receive such an insertion, otherwise 
the loop will give factors of m^. Simple power 
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counting shows that, for degenerate quarks, the 
diagrams which remain are those with any num- 
ber of rl loops coming out of a single mass-term 
vertex, as illustrated by Figs, lb and c. These 
are easily summed. Second, other diagrams in- 
volving more than one /x vertex can be shown to 
lead to non-leading logarithms using Weinberg's 
theorem 

(a) (b) (c) 



Figure 1. Contributions to tyIt^. The cross repre- 
sents the ry'^ vertex, the circle the \x term. 



I wish to recast the diagram summation into 
a RG framework. Start at a scale Ag « IGeV 
at which the coupling is /i(Ao) = /xq. Reduce A, 
varying /i(A) so as not to change amplitudes with 
small external momenta. When A drops to m^' , 
the loop diagrams cease to generate logarithms, 
and ^(A) no longer "runs". At this point use the 
tree level formula to determine m,,' : 

= m\, = 2[i(m^-)mq . (6) 

A crucial check on this program is that the RG 
equation does not depend on the choice of phys- 
ical amplitude used to determine it. The two- 
point, four-point, six-point, etc. amplitudes all 
yield the same result. This ensures that £bg re- 
tains its chiral invariant form. 

To obtain the RG equation, consider the two 
point amplitude. As already noted, the diagrams 
which contribute at leading order are those of the 
type shown in Fig. 1, but with any number of ry' 
loops. Let An be the contribution of the diagram 
with n 77'^ vertices. Each loop gives a factor of 
(51n(A^/m^,), and there is a combinatoric factor 
of 1/n!, so 

A„ = 2M(A)m,[(51n(AVm2,)]"/n! . 

Thus if we change In(A^) by (iln(A^), 

dAn = (5yl„_idln(A^) -I- A^dXwy. ; A_i . 



Summing over n we see that the total amplitude 
will be unchanged if /x(A) changes according to 

dln^ = -5dln(A2) ^ ^(A) = [iQ^t^H K^)^ . (7) 

The self-consistent equation, Eq. ^ becomes 

ml = ml, 2^iom,iAl/ml,y . (8) 

The solution is Eq. ||, except that we now know 
that /i is to be evaluated at the scale A. 

More generally, all amplitudes coming from the 
mass term in Cbg, ^^nd involving momenta less 
than m^/ , will come with an overall factor of 

^iim^,) = M^l/m^,)' . (9) 

Eq. ^follows from this result. 

This derivation is unusual in that one must 
work to all orders to obtain the RG equation. No 
extra diagrams are summed up when solving the 
equation. Furthermore, the external momentum 
does not act as an infrared cut-off on the diagram 
since it does not flow through the 77' loops. Never- 
theless, I think it is useful to sec the result emerge 
in this way. For one thing it may simplify the 
generalization to non-degenerate quarks. It also 
emphasizes that the non-analytic dependence on 
m'^ comes from the infrared part of the integrals. 

In Ref. I suggested removing the non- 

analyticities by assuming a non-analytic relation 
between the quark mass in the quenched chi- 
ral Lagrangian and that in quenched QCD. I no 
longer think that this is an option. The two La- 
grangians are matched at A ~ IGeV, in such a 
way that their infrared behaviors agree. Thus 
infrared divergences leading to non-analyticities 
cannot enter into the relationship between quark 
masses. 

3. COMPARISON TO DATA 

Present quenched data show no evidence for the 
divergences of Eqs. and ^. This is illustrated 
with the data of Ref. g in Fig. 2. Eq. | with S = 
0.2 gives the dashed line. (A is chosen to match 
the data at w 0.03.) One is free, however, 
to add terms linear and quadratic in m^. With 
coefficients of typical size, it is possible to make 
the theoretical curve quite flat, as shown by the 
solid line. Thus this data set probably cannot 
rule out S « 0.2. However, new results on the 
spectrum presented at this conference are likely 
to give a much smaller bound on S. 
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Lagrangian this would require a large value of the 
coefficient A mentioned above. 

If S is very small, then the problems with 
the quenched approximation are pushed to small 
masses. This may mean that they can be ignored 
for practical purposes. A direct calculation of 6 
would be very useful; another method for doing 
this is to study finite volume effects (^Q. 
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Figure 2. \n{m'^/mq) vs. (in lattice units). 



4. CONCLUSIONS 

I can think of three explanations for the lack of 
numerical evidence for the divergent terms. (1) 
The effect could be present (with S w 0.2) but 
hidden by a conspiracy of higher order terms. 
This I think is unlikely given the latest data, but 
I await with interest the results. (2) The theoret- 
ical analysis could be wrong. This seems unlikely, 
for the essential physics is in the 1-loop calcula- 
tion, and this requires only the presence of the 
light rj' and of the rj'^ vertex. See also the talk 
of Golterman (3) The effect is present but 
(5 <C 0.2. As pointed out to me by Fukugita, Itoh 
et al. attempted to directly measure ttiq with Wil- 
son fermions |^ . They found, with poor statistics, 
a number ~ 6 times smaller than (900MeV)^. 

I presently favor the third option. The only 
doubt I have stems from the following fact: to 
calculate mg one looks at the disconnected part 
of the pseudoscalar two point function, and finds 
the residue of the double pole at p"^ = — m^. 
The same two point function at p = is propor- 
tional to the topological susceptibility, x If 
the double pole gives the dominant contribution 
to the correlator at p = (with no momentum 
dependence in the residue), one can show that 
f^rriQ = 6x- Present numerical values of x then 
imply Too ~ 900 MeV. For toq to be much smaller 
than this the residue of the pole must have con- 
siderable momentum dependence. In the effective 
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